Quantum rearrangement and self-consistent BCS-BEC crossover thermodynamics 
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Based on previous works, analytical calculational procedures for dealing with the strongly inter- 
acting fermions ground state are further developed through a medium dependent potential in terms 
of the Bethe-Peierls contact interaction model. The methods are exact in the unitary limit regime 
and they lead to the self-consistent equations analogous to the Hartree ones. The single particle 
energy spectrum rearrangement effects on the thermodynamics due to the Hugenholtz-van Hove 
theorem constraint are addressed. These effects lead to an additional instantaneous correlation po- 
tential contribution to the system physical chemical potential and pressure, i.e., equation of state, 
which enforces the classical thermodynamic consistency. The Dyson-Schwinger equations represent 
implicitly the various Bethe-Goldstone expansion ones. In a thermodynamically self-consistent way, 
the universal dimensionless factor is analytically calculated to be ^ = |, which defines the ratio of 
the unitary fermions energy density to that of the ideal non-interacting ones at T = 0. 
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I. INTRODUCTION 

The strongly interacting matter physics serves cur- 
rently as a pivot topic in the fundamental many-body 
theory or statistical mechanics context. An interesting 
observable that characterizes this kind of physics is the 
unitary limit thermodynamics. 

Taking the Feshbach resonance techniques, it is now 
possible to magnetically tune the two-body inter-particle 
interaction strength. Increasing the scattering length of 
fermions from —oo to -|-oo resulting in bound boson sys- 
tems can exhibit the crossover from Bardeen-Cooper- 
SchriefFer(BCS) to Bose-Einstein condensation(BEC). 
This will lead to the universal thermodynamics of the 
Bethe-Peierls contact interaction fermions system with 
a zero-range, infinite scattering length interaction[l-ll]. 
To solve the universal resonant superfluidity thermody- 
namic problem, at least in principle, more theoretical ef- 
forts are urgent to understand the detailed dynamics by 
going beyond the mean field theory [12]. 

At the crossover point, the dimensionless coefficient 
^ relates the total energy E/N = Cf^/ with the 
Fermi kinetic energy e/ = k'j/{2m). Here, m is the 
bare fermion mass while fc^ is the Fermi momentum. 
The determination of the fundamental constant ^ as 
a challenging many-body topic attracts much attention 
theoretically/experimentally[8, 13-16]. It is presently ac- 
cepted to be in the region 0.40 — 0.46. Compared with 
the old values, the updating experiments converge more 
towards C ^ 0.42 - 0.46[17-19]. 

The ultra-cold fermion atoms gas physics offers a plau- 
sible perspective in looking for the general nonperturba- 
tive statistical field theory methods. In the strongly in- 
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teracting regime kf\a\ » 1, especially right at the uni- 
tary limit point a ~ ±oo, the infrared singular nature of 
the two-body scattering amplitude /o(/c,a) = i/k (with 
k being the relative wave-vector magnitude of the col- 
liding particles) rules out any conventional perturbative 
expansion techniques. It is generally believed that there 
is not a simple analytical method to compute the uni- 
versal coefficient ^, although this topic stimulates many 
growing arduous theoretical efforts. By noting that an 
exact many-body theory concerns the behavior at uni- 
tarity is presently not available, it is not surprising that 
the existed analytical results about the universal coeffi- 
cient ^ differ from each other considerably ~ 0.2 — 0.6). 
It is worthy noting that there are still uncertain discrep- 
ancies among the powerful quantum Monte-Carlo works 
themselves as well as the experimental ones existed in 
the literature[8, 14, 15] [17-19]. 

In Refs.[20], we make an almost first principle de- 
tour to attack the universal constant ^ with the gen- 
eral Dirac quantum many-body formalism. The analyti- 
cal circumvention is motivated by the strongly instanta- 
neous Coulomb correlation thermodynamics discussions 
in a compact nuclear confinement environment [2 1-23] 
based on the relativistic continuum nuclear many-body 
framework[24, 25]. The infrared correlation character- 
istic in the unitary limit promotes us to use the uni- 
versal strongly instantaneous long range interaction to 
attack the infinite contact interaction thermodynamics. 
The obtained analytical result ^ is comparable with other 
theoretical quantum Monte-Carlo attempts. It is derived 
through a medium dependent interaction by accounting 
for the interior correlation with a counteracting or anti- 
correlation external source static potential. 

However, the auxiliary quantities pressure P and 
chemical potential // do not obey the profound universal 
hypothesis in the unitary limit with T = 0[4]. According 
to the universal hypothesis at unitarity with T = 0, the 
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ratio of pressure to energy density should be still that 
of non-interacting fermion gas. Therefore, the physical 
pressure should be P = ^Pideal instead of P = ■^Pideal ob- 
tained in Refs.[20] and consequently affects the important 
sound speed. The chemical potential should be ^ = ^Cf 
instead oi p. = /x* +27rae//n/m = e f — n^n/ {mkf) = ie/. 
Intriguingly, the physical Helmholtz free energy density is 
not affected. At T = 0, it can be / = —P + jln = ^^nej. 
Therefore the underlying encrgy(-momentum) conserva- 
tion law is still exactly ensured even with the respective 
discrepancy of the auxiliary pressure and chemical poten- 
tial. This observation instructs us to reanalyze the phys- 
ical pressure and system chemical potential carefully. 

The Hugenholtz-van Hove(HvH) theorem characterizes 
the fundamental thermodynamic relation between energy 
density and pressure[26]. It is well known that any ef- 
fective medium dependent interaction with Brueckner- 
Bcthe-Goldstonc(BBG) techniques can easily violate this 
relation. The standard many-body technique to remove 
this tortuous difficulty is considering the quantum rear- 
rangement effects of the single particle energy spectrum 
on the physical chemical potential and prcssurc[27 30]. 
The role of the single particle energy spectrum shift will 
be explored in the strongly interacting unitary limit. 

To elucidate the physical chemical potential and pres- 
sure(equation of state) at both finite temperature and 
density transparently with a more concise non-relativistic 
recipe for the ultra-cold atomic physics is the present 
main motivation. Up to now, quantum statistical field 
theory of the strongly interacting matter at finite tem- 
perature and density is still a challenge. More concretely, 
the exact determination of equation of state plays an im- 
portant role in looking for the detailed knowledge about 
the collective excitation modes measured with the ultra- 
cold atomic physics experiments. 

The present work is organized as follows. In Sec. II, 
the adopted low energy effective theory formalism and 
corresponding nonperturbative approach for the strongly 
interacting unitary gas thermodynamics are described. 
Through taking into account the single particle energy re- 
arrangement effect, the unitary fermions thermodynam- 
ics with the dynamically self-consistent Dyson-Schwinger 
equations is presented in Sec. III. Moreover, the univer- 
sal thermodynamics at unitarity with T = as well as the 
concrete comparisons with existed results is discussed in 
Sec. IV. Section V includes the discussion and prospec- 
tive remarks. The summary conclusion is made in the 
final Sec. VI. 

In this work, the calculations are performed in the non- 
relativistic contact interaction formalism. Throughout 
the paper, the natural units h = ks = ^ are employed. 



II. FORMALISM 
A. In-medium effective action 

The usual bare Bethe-Peierls Hamiltonian with a 
global U{1) symmetry addressing the many-body ground 
state energy can be found in the literature. The in- 
medium behavior associated with the many-body corre- 
lation characteristic is the key, while a nonperturbative 
approach is crucial for the strongly interacting unitary 
regime physics. Motivated by previous works, we work 
with the medium-scaling universal functional version 

f 

+^%^ J d^xr„ixW0ix)M^)Mx). (1) 

In Eq.(l), m is the bare fermion mass with a, (3 
, I representing the (hyperfine-)spin projection. The 
physical many-body renormalized interaction potential 
U*j: J [n] is a functional of density and temperature, from 
which one can define the in-medium S-wave scattering 
length aeff 

t/:>,[n,T].l^I^^. (2) 

The bare interaction matrix /7*^j[n, T] Uq is directly 
related with the free-space scattering length a with the 
low energy Born approximation for the zero range contact 
interaction. 

Considering the medium influence of the surrounding 
environment on the single particle energy spectrum prop- 
erties, we have included the rearrangement term oc /Ura[?^] 
in Eq.(l). It is introduced through the auxiliary quantity 
fla indicated by the tilde symbol fla [n] = fia— lira [n] ■ It 
means that the collective correlation effects arc further 
accounted for as a single particle potential formalism in 
the spirit of density functional theory [27, 31], while the 
true Lagrange multiplier is /ia- As a functional of den- 
sity and temperature, its concrete medium dependence 
can be derived analytically from the thermodynamic re- 
lation for the fixed particle number density. 

B. Effective potential determined by the 
counteracting environment 

Theoretically, the derivation of the many-body renor- 
malized effective interaction U*^^ is beyond the bare 
Bethe-Peierls contact interaction Hamiltonian itself, i.e., 
it should be derived from the more underlying physical 
law[27]. 

Firstly, from the point of view of the contact interac- 
tion fermions physics, it is extremely instructive to recall 
the primary intermediate vector boson (IVB) hypothesis 
in weak interaction theory. It is well known that the lo- 
cal intermediate vector boson theory is related with the 
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current-current (CC) contact interaction version through 
the corresponding connection: g^/rri^ = In the 

low energy limit k ^ 0, the two IVB and CC theories 

arc ingraincdly identical. To model the low energy uni- 
tary limit physics, we "let" g^, /rriyy — > oo by taking gw 
to be an arbitrary large charge or with an arbitrary small 
mass gap myy- Consequently, the counteracting environ- 
ment associated with an additional U(l) conserved "elec- 
tric charge" can be introduced to attack the challenging 
infinity. Constrained by the opposite charge neutraliz- 
ing bac;kground, the effective interaction is found to 
appear as an instantaneously anti-screening formalism. 
This particular analogy and/or discursion makes the cou- 
pling constant gw enter simultaneously in the denomina- 
tor and numerator of fractions in the relevant thermody- 
namic expressions [20]. The assumed electric charge gw 
is reduced by the physical constraint or charge neutral- 
izing stability condition and its concrete magnitude does 
not affect the final unitary physical results. Therefore, 
there are not any regularization or expansion parameters 
violating the intriguing unitarity in gauging the infrared 
singular divergence. This spirit will be further explored 
in this more concise non-rclativistic recipe. 

Moreover, corresponding to counteracting the devi- 
ation tendency of charge neutralizing, there exists a 
similar general Le Chatclier's principle in thermody- 
namics. This profound stability principle accounts for 
frustrating an instantaneous departure from equilibrium 
with an alternating minus function of the surround- 
ing environment [32]. First of all, as we will sec, the 
counteracting-antiscreening formalism makes it possible 
for us to use the algebra equations to characterize the 
classical thermodynamics with the significant quantum 
correlation effects, but of crucial importance is that 
the many-body renormalization effects are introduced 
through a mirror background with an alternating minus 
function to avoid the theoretical double counting trouble. 

The following calculations based on Eq.(l) can be char- 
acterized by the comprehensive Feynman diagrams pre- 
sentation with the generalized coupled Dyson-Schwinger 
equations(DSEs) formalism as indicated by Fig.l. 

The polarization tensor n*(a;„,p) is calculated 
through the non-relativistic random phase approxima- 
tion(RPA) with the full fermion propagator ^(a;„,k) as 
described by Fig. lb 

a;„,k-|-p) 

fk+p — /k 



V 2m 2m 



(3) 



with Ji^ = J (i'^k/(27r)'^. Here, /k is the Fermi-Dirac dis- 
tribution function given explicitly in Eq.(16). 

After employing the analytical continuation for the 
Matsubara frequency ojn according to iu>„ w + irj, the 
n*(a;,p) determines the dispersion relation of the exci- 
tation mode A (the composite boson/dimer) in- medium 
propagating. The polarization function obtained in this 
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FIG. 1: Generalized self-consistent Dyson-Scliwinger equa- 
tions incorporating the Landau pole contributions associated 
with the Fermi surface (anti-)correlation instability [20]: (a). 
The full fermion propagator determined with antiscreening in- 
teraction potential in instantaneous approximation U^f / [n] = 
Uo/{l + n*(0, IpI 0)Uo); (b). The boson propagating cal- 
culated with the full fermion propagator. 



way is the standard in-medium Lindhard density-density 
or spin-spin response function for the one component 
fermions gas [12]. In this work, we care about the low en- 
ergy long wavelength thermodynamics of the fully pair- 
ing correlation occasion, which is related to the static 
infrared limit 

n*(0,|p|^0) = -I^ = i(!^)3/^Liv.[.].(4) 

The mjj appearing in Eq.(4) is the physical non- 
relativistic Debye screening mass squared with unit elec- 
tric charge and reflects the quantum fluctuation contri- 
butions consequently. Its magnitude is the spin-spin or 
density- density correlation response susceptibility in the 
static infrared limit. At T 
the familiar formalism 



0, mjj can be reduced to 



2 I ^f"^ 



(5) 



Different from the multiplier-system physical chemical 
potential /i, the symbol A employed in above equations is 
the effective(total) chemical potential. It characterizes 
the additional many-body collective correlation effects 
and will be defined explicitly in the following discussion 
for the full fermion propagator. Here and afterwards, the 
function identifier Li^ [y] stands for the standard Fermi in- 
tegral mathematica polylogarithm function polyLog[i, y]. 
Meanwhile, we will use the shorting notation y = — eT 
for the sake of conciseness throughout the paper. 

With the polarization hmction H*(0, |p| 0), the re- 
arrangement interaction matrix is derived to be 



U:ff[n,T] 



Uo 



i + n*(o,|p|^o)c/o 1 



-Uo 



(6) 
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At T = 0, it can be reduced to 



U. 



effl 



l-N{ef)Uo' 



(7) 



The A^(e/) = kfm/{2TT'^) is the un-pcrturbatcd den- 
sity of states on the Fermi surface for one component 
fermions[12], which is one half of the physical Debye 
screening mass squared. Without a loss of generality and 
conformity, the screening mass squared notion rnjj in- 
stead of the density of states will be used to characterize 
the general finite temperature thermodynamics. 

Compared with the conventional screened interaction 
customarily applied in the weak coupling perturbative 
loop diagram ring/ladder resummations 



C^md(w„,p) 



f/n 



i-n*K,p);7o' 



(8) 



there is a key dynamically mirrored minus sign difference 
in the denominators of Eq.(6) and Eq.(8), which leads to 
the quite different physical motivations and calculational 
details. 



III. THERMODYNAMICS WITH 
GENERALIZED DYSON-SCHWINGER 
EQUATIONS 

A. Full fermion propagator and energy density 

The full fermion propagator ^(a;„,k) given with the 
instantaneously anti-screening interaction matrix tells us 

all the information of the ground state energy or quasi- 
particlc distribution function. It is calculated with the 
proper self-energy E* through the instantaneous approx- 
imation as indicated by Fig. la 



Gt = Got + ^OT^*^T' 



(9) 



n 



S*(0,0) = U:ff{0, |k| ^ 0)ni = Kffin]-. (10) 

Here, n is the total particle number density n — n-^ +ni = 
2nij for the unpolarized/fully pairing symmetric trapped 
dilute atomic Fermi gas. 

The self-energy S*(0, 0) is four-momentum indepen- 
dent in the instantaneous approximation. It means that 
the self-energy will shift the chemical potential of the par- 
ticle distribution function, i.e., the correlation contribu- 
tions to thermodynamics are characterized through the 
phase space energy-momentum distribution deformation. 

From the full fermion propagator with Eq.(l), the par- 
ticle number and energy densities are given by [33] 



N 
V 

E 
V 



T I ^e-"''tra(c^„,k), (11) 
^/E^""''^(^^" + 1^ +/i)tra(a;„,k).(12) 



Although wc once again consider the linear contact 
formalism, the lowest order diagrams characterize the 
thermodynamics with the highly nonlinear /turbulent 
contributions of the density and temperature fluctua- 
tions due to taking the medium dependent interaction 
potential[27]. 

After summing the Matsubara discrete frequency a;„ 
with the residual theorem, the particle density n = N/V 
and energy density e = E/V are reduced to the compact 
formalisms 



n 



= 2^/k = -2(^)3/2Li3/2[y], (13a) 



e = 2 



H n{l,X) 

2m m 



= -3T(!^)3/^Lv,M+ ^ 
with the entropy density s = S/V 



^ff 2 



(13b) 



= -2 / [/k In /k + (1 - /k) ln(l - /k)] 

Jk 

= U^)'^' (2ALi3/2[t/] - 5TLi5/2[t/]) . (14) 



T' 27r 

The in-medium effective scattering length is 

a 



^eff 



1 27rm?, 
1 7:r^a 



(15) 



In deriving above equations, we have defined the quasi- 
particle Fermi-Dirac distribution function with the to- 
tal/effective chemical potential A through Fig. la with 
Eq.(9) 



/k 



k2 



■^) + l 



A =/i- S*(0,0), 



(16) 
(17) 



with P = 1/T being the inverse temperature. 

The expressions for the entropy density Eq.(14) and 
the first energy density term in Eq.(13b) are similar to 
those of the non-interacting fermion gas. The total chem- 
ical potential A is determined for the given particle num- 
ber density n and temperature T. The above algebra 
equations are completely closed. 



B. Physical pressure and system chemical potential 
with recU-rangement 

From the Helmholtz free energy density F/V = f = 
e — Ts, the pressure is given by 



dv 



T,N 




di 



-(18) 



dn 



At T = with entropy density s = 0, P is reduced to the 
important formalism 



P = n' 



d_ 

dn Vn 



(19) 
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The equation Eq.(19) determines the fundamental rela- 
tion between the energy density and pressure according 
to the HvH theorem[26] , which plays a vital role for char- 
acterizing the symmetric nuclear matter saturation prop- 
erty. 

As widely discussed in the nuclear many-body liter- 
ature, the HvH theorem constraint will lead to the re- 
arrangement terms in the pressure and chemical po- 
tential dealt with any kind of in-medium interaction 
potential[27-30]. Presently, their exact analytical expres- 
sions can be derived in a thermodynamic way instead of 
doing the numerical interpolation[34]. 

Let us further address the pressure P. From Eq.(18) 
given with the medium dependent interaction potential, 
one can have 



m 



n^ + Pideai{T, A), 



(20) 



with 



Pideal {T, A) 



fmT\ 



3/2 



Li5/2M. (21) 



In Eq.(20), the first term represents the expected direct 
contribution of the induced interaction, while the last 
one results from the Fermi kinetic energy. The essential 
second term is due to the rearrangement contributions, 
which can be easily neglected. It is worthy noting that 
this rearrangement term is vanishing within the standard 
mean field theory approach [34]. 

The discussion for the physical system chemical poten- 
tial/multiplier is in line with that of pressure 



dF 
dN 



T,V 




d{e - Ts) 
dn 



■ (22) 



Its anal3^ical expression can be further reduced to 
M = + C{T, A) 



+ A. (23) 



There is an additional(second) term in the system chemi- 
cal potential resulting from the intriguing rearrangement 
effect of the single-particle degrees of freedom analogous 
to that of pressure. 

Explicitly, the rearrangement terms in the pressure 
Eq.(20) and chemical potential Eq.(23) are similar. For 
the general finite temperature occasion, the rearrange- 
ment factor C(T, A) appearing in Eq.(20) and Eq.(23) is 
given by 



C(T, A) = ruD 



dn 



Li-i/2M 
2TLii/2[i/]' 



(24) 



From the pressure and chemical potential analytical 
expressions, one can see that the rearrangement effects 



play the positive "repulsive" role. These effects pre- 
vent the system from collapsing in the strongly coupling 
attractive unitary limit. Furthermore, the in-medium 
spontaneously generated rearrangement effects are ex- 
actly counteracted (easily neglected) in the Helmholtz 
free energy density / = e — Ts = —p + fj,n. This re- 
flects an important physical fact that the in-medium off- 
shell(spontaneous Lorentz violation) formalism can ap- 
proach the quantum correlation effects on the univer- 
sal thermodynamics, while the thermodynamics first law 
constraint ensures the energy conservation. Exactly, the 
classical thermodynamic consistencies are satisfied with 
the generalized self-consistent DSEs. 



IV. UNIVERSAL THERMODYNAMICS IN THE 
UNITARY LIMIT WITH T = 

Prom the general analytical expressions, one can give 
the T = results with the in-medium interaction 



Uo 



eff 



kfm 1 



"'eff 



1 



2fefa 



With n = fc^/(37r^) and the rearrangement factor 
C{ef)=mn^ = — 



(25) 



(26) 



the physical internal energy density Eq.(13b), pressure 
Eq.(20) and the system chemical potential Eq.(23) are 



E 


1 


2akf 


3 

+ 5- 




Ncf 


37r 


— 2akf 




PV 


1 


2akf 




' 2akf 




37r 


— 2akf 




^TT — 2akf 


a. 


2 


2akf 




' 2akf 




37r 


— 2akf 


^TT — 2akf 



2 

+ 5' 
+ 1. 



(27a) 
(27b) 
(27c) 



In the strongly interacting universal unitary limit 
\a\ oo, the ratios of energy density, pressure and chem- 
ical potential to those of non-interacting fermion gas are 
^ = |. This ^ = I is exactly consistent with the existed 
theoretical rcsults[8, 13, 15] and in agreement with the 
updating experimental ones[17-19]. In the weak coupling 
limit, the analytical thermodynamic expressions read- 
ily return to the lowest order mean field theory results 
which neglect the quantum fluctuation/correlaion contri- 
butions. 

At unitarity and with T = 0, the pressure is 



P = 



45m 



{3n')y\^^'^' = ^PUeal, 



with which the sound speed is calculated to be 



dP 2 r- 

man 3 



(28) 



(29) 
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The vpG = Y 1^/ is that for the ideal non-interacting 
fermion gas with vj = kf/m, being the Fermi speed. A 
detailed information of the sound speed of the homoge- 
neous system is very crucial for determining the frequen- 
cies of collective modes in the realistic trapped systems. 

With T = 0, the ratio of pressure to binding energy 
is still I of the non-interacting fermions gas in the uni- 
tary limit as pointed out in the literature [4]. The HvH 
theorem ensures the thermodynamic consistency dealt 
with the in-medium interaction potential in addressing 
the strongly interacting physics. The same ratio is coun- 
terintuitive to some extent, which indicates that the uni- 
versal thermodynamic properties with the infrared long 
range quantum fluctuation/correlation in the strongly in- 
teracting limit are extremely profound. 

The numerical results for T ^ can be readily given 
in the unitary limit regime with the explicit analyti- 
cal expressions, from which one can see that the collec- 
tive strongly correlating and rearrangement effects influ- 
ence the chemical potential and pressure simultaneously. 
However, the internal energy and Helmholtz free energy 
densities do not directly rely on the intriguing rearrange- 
ment effect even with a careful numerical check. 



V. DISCUSSION AND PROSPECTIVE 

In addressing the unitary fermions physics, the essen- 
tial task is how to incorporate the significant collective 
density and temperature fluctuation/correlation contri- 
butions into the thermodynamics in a reasonably con- 
trollable way. 

The physical connotations are crystal clear in this ana- 
lytical attempt: On one hand, the many-body surround- 
ing environment modulates the bare two-body interac- 
tion, which is characterized by the in-medium renormal- 
ized interaction matrix U*jj[n,T]. On the other hand, 
the collective correlation effects further shift the single 
particle energy spectrum, which is described by the cor- 
responding single particle correlation potential Hr[n,T]. 
The crucial energy-momentum conservation law is guar- 
anteed by the simultaneous consideration of these two 
collective correlation effects. The fundamental thermo- 
dynamic consistency is the basic validity requirement of 
any scientific theoretical approaches, which can be easily 
neglected. 

The generic infrared singularity of the bare two-body 
scattering amplitude at unitarity provides the scaling ba- 
sis, which leads to the thermodynamic universal hypoth- 
esis. This scaling property implies that the thermody- 
namic expressions can be very simple in principle, i.e, 
the analytical formulae should be similar to those for the 



ideal non-interacting fermions gas. In turn, this infrared 
singularity makes any potential theoretical approaches 
nontrivial. Essentially, the unitary characteristic rules 
out any naive expansions or conventional perturbative 
loop diagram ring/ladder resummation calculations. 

In this work, we have obtained the T = dimension- 
less universal coefficient ^ = | for the D = 3 unitary 
fermion gas. It is in the regime of the reasonable simu- 
lation/experimental values 0.42 — 0.46. However, as dis- 
cussed recently in the literature[35, 36], it is yet far from 
the ending of looking for its realistic exact value. 

The exact determination of this fundamental constant 
will influence the detailed knowledge about the collective 
excitation modes measured in the realistic ultra-cold uni- 
tary atomic fermions. Consequently, more and more the- 
oretical/experimental efforts are in fact most expected to 
make it certain precisely. To understand the comprehen- 
sive thermodynamic properties exhibited by these novel 
trapped strongly interacting quantum systems, there is 
still a faraway route before us to constitute a consistent 
microscopic theory and corresponding analytical calcula- 
tional framework. 



VI. CONCLUSION 

In summary, based on the Dyson-Scwhinger equations 
formalism, the ground state of the D = 3 ultra-cold uni- 
tary fermion gas has been calibrated in detail in terms of 
the quantum Ising-like Bethe-Peierls contact interaction 
model. The universal factor is calculated to be ^ = |. 

In addressing the universal resonant thermodynamics, 
the medium dependent interaction through a dynami- 
cally twisted minus sign caused by an assumed mirror 
background makes it possible for us to incorporate the 
Landau pole contributions with the Fermi surface insta- 
bility. With the spontaneously generated rearrangements 
due to the single particle energy spectrum shift, the ad- 
ditional instantaneous correlation potential ensures the 
classical thermodynamic consistency. Finally, we want 
to emphasize that there are not any infrared and/or ul- 
traviolet cutoffs or regularization and/or expansion pa- 
rameters obviously violating the significant unitarity, i.e., 
the bare Uo drops out in the unitary limit. 
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